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Abstract 

0\ ■ 

>: 

Oh' We have estimated the dispersion of the inferred relative frequencies of Type- la and Type- II supernovae 

o ■ 

(Ni^/Nii) in our Galaxy by fitting the numerical results of supernova nucleosynthesis to the solar-system 
abundances. The ratio Ni a /Nn is estimated to be 0.056-0.14, which is consistent with the observation, if 

> : 

the model of Woosley and Weaver (1995, WW95) is adopted for Type-II supernovae (SNe II). On the other 
hand, the upper limit of Ni a /Nu becomes too large for the model of Hashimoto (1995, Ha95). However, 
Ha95 can fit the solar values better than WW95 as far as the abundant nuclei are concerned. These results 
mean that Ha95 can reproduce the main nuclei of the solar-system well and that WW95 can reproduce 
the solar values over a wide mass-number range. We also note that Ni a /Nu tends to become smaller if the 
delayed detonation model is adopted for Type-la supernovae (SNe la). The dispersion of Ni a /Njj obtained 
in the present investigation is larger than that concluded by Tsujimoto et al. (1995), which will have an 
influence on the estimate of the average life time of SNe la's progenitors and/or the star- formation rate 
history in our galaxy. 
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Supcrnovae play an important role in ejecting heavy elements of A > 12 produced during stellar evolution and/or 
an explosion. There are two types of supernovae: Type-I and Type-II supernovae. Type-la supernovae (SNe la) are 
thought to be thermonuclear explosions of accreting white dwarfs (e.g., Nomoto et al. 1994). On the other hand, 
Type-II supernovae (SNe II) and Type-Ib/Ic supernovae (SNe Ib/Ic) are thought to be collapse-driven supernovae 
(e.g., Bethe 1990). Many calculations of them have been performed and their chemical compositions have been 
calculated. 

As for SNe la, the defraglation model has been thought to be more favored than the detonation model, since 
56 Ni is over-produced in the detonation model. However, concerning the amount of 56 Ni, the defraglation model 
can explain the spectrum of SNe la fairly well (e.g., Nomoto ct al. 1997a). In particular, a model W7 (Nomoto et 
al. 1984) has been considered to be 'standard', and the calculated chemical composition has been used widely to 
explain the solar-system abundances. As for SNe II, two groups have systematically calculated the nucleosynthesis 
(e.g., Woosley, Weaver 1995; Hashimoto 1995). Although their calculations assumed different physical process such 
as the criterion of convection, the 12 C(a,7) 16 rate, which is crucial for constructing stellar models, and initial 
shock-wave conditions, which determines the peak temperature during the explosion, they succeeded in explaining 
the chemical compositions of SNe II with a fair degree of precision. In fact, the solar-system abundances are well 
reproduced within a factor of 2-3 by the combination of their results for a wide range of mass numbers, except for 
some nuclei (Hashimoto 1995; Timmes et al. 1995.) 

By using the results of W7 and Ha95, the relative frequencies of SNe la and SNe II in our galaxy was obtained to 
be Aia/iVn = 0.15 (Tsujimoto et al. 1995), which is consistent with the fact that the observed estimate of the SNe 
la frequency is about 10% of the total supernovae occurrence (van den Bergh, Tammann 1991). However, there is 
one question in their calculations. They concluded that the contribution of SNe la to the solar-system abundances 
is in the range r — 0.09 ± 0.01 against possible uncertainties in their analysis (see subsection 2.2). However, they 
arbitrarily selected 31 species to be fitted, and neglected other nuclei when they derived the proper value of r. We 
think that it is worth examining whether the dispersion of the Ni^/Nu changes or not by the choice of the nuclei 
to be fitted. Moreover, we should examine whether the proper value of r can be really predicted from numerical 
calculations without such an arbitrariness. 
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Another motivation for the present study is the model dependence of the prediction of Ni a /Nu. Recently, based 

on the progress of numerical calculations for supernovae, the chemical compositions of their ejecta have been recal- 
culated. In particular, the effect of the delayed detonation has been taken into consideration for SNe la. As for the 
W7 model, the flame speed is assumed to be relatively high, which is inconsistent with that of multi-dimensional 
calculations; the simple defraglation can not generate an energetic explosion, or can not lead to an explosion in 
the multi-dimensional calculations (Arnett, Livne 1994; Khokhlov 1995). To solve this problem, a mechanism of 
delayed detonation is proposed by Khokhlov (1991). Recently, nucleosynthesis in SNe la was recalculated based on 
the assumption of delayed detonation (Nomoto et al. 1997b), and has presented a more excellent explanation for 
the observation of SNe la. For example, the velocity profiles of intermediate nuclei, such as Si and S, become wider 
compared with that of the W7 model, which is consistent with the observations. We also note that the chemical 
compositions of SNe II are different between Ha95 and WW95 (Hashimoto 1995; Woosley, Weaver 1995). 

To investigate the capability of reproducing the solar-system abundances with using the avairable results, we will 
explore the model dependence and derive the dispersion based on numerical calculations. It is important to estimate 
the dispersion, since the change in the relative frequency of SNe la and SNe II will have an influence on the estimate 
of the average lifetime of SNe la's progenitors and/or the star- formation history in our galaxy. 

We describe our calculation method in section 2. The results are presented in section 3. A summary and a 
discussion are given in section 4. 

2. Models and Method of Calculations 

2.1. Models of SNe la and SNe II 

As for the composition of SNe la, we use the result of W7 for the simple deflagration model and WDD2 for 
the delayed detonation model, respectively (Nomoto et al. 1984; Nomoto et al. 1997b). We also use the results 
of Ha95 and WW95 for SNe II. More precisely, we use for Ha95 the results of helium-star models with masses of 
M a = 3.3, 4, 5, 6, 8, 16, and 32, which approximately correspond to main-sequence stars of M ms = 13, 15, 18, 20, 25, 
40, and 70 M Q , and for WW95 S13A, S15A, S18A, S20A, S25A, and S40A models. The averaged synthesized yields 
in 10-50M Q stars are obtained with a weight of Salpeter's initial mass function, as follows: 



4 



S. Nagataki, M. Hashimoto, and K. Sato 



[Vol. 



The differences in the calculated chemical compositions are shown in figure 1. The open circles denote W7/WDD2 
and the dots represent WW95/Ha95. The abscissa represents the mass number. The plotted 46 nuclei, which are 
the same nuclei as those in table 1 of Tsujimoto et al. (1995), are summarized in table 1. 

2.2. Calculation Method for the Solar- System Abundances 

In an analysis of the reproduction of the solar-system abundances, we adopt the same method with Tsujimoto ct 
al. (1995). At first, as the sum of Mij a we define M\ a which are the total heavy element mass of SNe la. We also 
define Mu in the same way. 

Next, we define the abundance pattern Xi clS 



which is to be compared with the solar Xi(&), defined as Z i /^2 i Zi, where Zi is the observed abundance of the i-th 
element per unit mass (Anders, Grevesse 1989). The most probable value of r — r p is determined by minimizing the 
following function (Yanagida ct al. 1990): 



where i runs over the heavy elements and their isotopes chosen in the minimization procedure. We summarize those 
elements in table 1. We note that r has a different meaning from that of the relative frequency Nj^/Nu- The relation 
between them is 

-/Via = u>u M u r p 

N u w Ia M Ia (l-r p )' U 
where wi a and Wn represent the mass fraction of heavy elements ejected into the interstellar gas from SNe la and 
SNe II, respectively. These values were estimated to be 0.27 and 0.22 in the solar neighborhood from the numerical 
calculation (Tsujimoto et al. 1995). We calculated g(r) in five cases. We used all 46 nuclei (Case 1). Next, we used 
the selected 31 nuclei (Case 2) and the most abundant 6 nuclei (Case 3). Then, 20 elements, each of which is a 
summation of the isotopic abundances, were examined (Case 4). Finally, the selected 14 elements were fitted (Case 
5). Furthermore, we also performed the same analysis for LMC and SMC with the same elements with Tsujimoto 



Xi = rM iM /M la + (1 - r)M i ,„/M„ (0 < r < 1), 



(2) 



n 




(3) 
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et al. (1995) to check our conclusion. The value of Mu/M\ a for each case is summarized in table 2. We note that 

Tsujimoto ct al. (1995) concluded that r p = 0.09 ± 0.01 and -/V Ia/ //V n = 0.15 in their analysis by adopting Case 2 

and Case 5 for the results of W7+Ha95. 

3. Results 

3.1. The Solar System Abundances Pattern 

We show in figure 2 the form of g{r) for four models: W7+Ha95, W7+WW95, WDD2+Ha95, and WDD2+WW95. 
The thick solid, long-dashed, short-long-dashed, short- dashed, and dot lines correspond to Case 1, 2, 3, 4, and 5, 
respectively. 

Generally speaking, g(r) has a minimum value at r ~ 0.1 in most cases. In particular, g(r) has minimum values 
for all cases if WW95 is adopted. On the other hand, g(r) is not sensitive to r for the W7+Ha95 and WDD2+Ha95 
models in Case 1 and 4. However, g(r) has a deep minimum value in the case of 3 if Ha95 is adopted. These results 
mean that WW95 can reproduce the solar values well in many nuclei and that Ha95 can reproduce them well for 
selected abundant nuclei. The values of r p for all cases are summarized in table 3. We should note that the range 
of r p is wider than that of Tsujimoto et al. (1995). 

We can see the effect of delayed detonation. The values of r p tend to become smaller if WDD2 is adopted. We 
also feel that the minimum values of g(r p ) seem to become smaller in WDD2. This means that the solar-system 
abundances are fitted better by WDD2 than W7. We examine this tendency by adopting our analysis to SMC and 
LMC in the next subsection. 

3.2. The SMC and LMC Abundance Patterns 

In subsection 3.1, we found that both r p and g{r p ) become smaller if WDD2 is adopted. We will check this 
tendency by adopting our method of analysis to other galaxies: SMC and LMC. As was done by Tsujimoto et al. 
(1995), 12 and 10 elements were taken into consideration for SMC and LMC, respectively (see table 1). The results 
are given in table 4. We found the same tendency of r p and g(r p ), which supports our conclusion. We also note the 
r p s obtained in this analysis tend to be larger than those of the solar-system abundances, which is consistent with 
Tsujimoto et al. (1995), even if WDD2 is adopted. 
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3.3. Relative frequency of SNe Ia/SNe II 

The ratio iV Ia /iVn depends on the three values of Wn/wiai Mu/Mi a , and r p . In our analysis, Wn/wi a is constant 
(subsection 2.2). WW95 tend to make M\\jM\ & smaller than Ha95. WDD2 also make it lower, though its degree is 
weak. The range of N\ a /Nu in table 3 is estimated from r p and M\\JM\ a for Case 1-5. The values in the parenthesis 
in table 3 are not so accurate, since g(r) is insensitive to r. We note that Ni a /Nu is estimated to be 0.056-0.14, 
which would be consistent with the observation if the model of WW95 is adopted. On the other hand, the upper 
limit of Ni a /Nu becomes too large for Ha95. We also note that Ni a /Nu tends to become smaller if WDD2 is adopted. 
Generally speaking, the range of Ni a /Nu is fairly wider compared with the conclusion of Tsujimoto et al. (1995). 

3.4-. Other Analysis 

The method described in subsection 2.2 is not the only one that can be used to determine r p . In this subsection 
we examine another method to obtain r p : \ 2 fitting. We use two elements, O and Fe, for the x 2 fitting. The r p is 
determined to be 7.5 x 10~ 2 -1.0 x 10" 1 , (6.8-9.1) x 10~ 2 , (1.1-1.3) x 10" 1 , and 9.8 x 10~ 2 -1.2 x 10" 1 (99% C.L.) for 
W7+Ha95, WDD2+Ha95, W7+WW95, and WDD2+WW95, respectively. However, we must say that all models 
are excluded unless proper nuclei are chosen for the fitting, since the uncertainty of the observation is much smaller 
than that of the theoretical prediction. That is why these ranges mentioned above are less persuasive than those 
discussed in subsection 3.1. In other words, we must wait for progress concerning theoretical calculations to use this 
method effectively. 

3.5. Present Situation of Numerical Simulations 

In figure 3, we show the solar abundance pattern Xi/xi(Q) (see, equation (2)) predicted by numerical calculations 
using r p of Case 2. The open circles represent WW95 and the dots denote Ha95. There are some nuclei which 
are less produced in the range of mass numbers 35-50. For example, Ha95 has a less-production problem of some 
elements, such as Ar, K, Ti. WW95 also has a problem of less production of 41 K, 44 Ca, and Ti. Cu and Zn are 
also less produced by Ha95. On the other hand, some nuclei in the range 50-65 are over produced. These problems 
remaine to be solved to predict r p more precisely. 
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4. Summary and Discussion 
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We have explored the dispersion of Ni a /Nu using the four models Ha95, WW95, W7, and WDD2. As for r p , 
g(r) clearly has minimum values if WW95 is adopted. On the other hand, g(r) is insensitive to r for the W7+Ha95 
and WDD2+Ha95 models in Case 1 and 4. However, g(r) has a deep minimum value in the case of 3 if Ha95 is 
adopted. These results mean that WW95 reproduces the solar values well in many nuclei and that Ha95 reproduces 
them well especially with regard to abundant nuclei. Among the most abundant 6 nuclei, the abundance of 56 Fe 
is different most between WW95 and Ha95. This would be because only WW95 determines the mass cut by a 
hydrodynamical calculation. This suggests the difficulty of the determining the mass cut from the hydrodynamics, 
even if the spherical explosion is assumed. Mu/M\ a tends to become lower for WW95 than for Ha95. Although 
WDD2 also make it lower, its degree is weak. 

As a result, the range of N\ a /Nu becomes 0.056-0.14 for WW95+SNe la , which is consistent with the observation, 
and 0.077-(0.61) for Ha95+SNe la. We also note that Nj^/Nu tends to become smaller if the delayed detonation 
model is adopted for Type-la supernovae. Since the delayed-detonation model may be a more realistic model than 
the simple dcfraglation model, this tendency should be stressed. Generally speaking, the dispersion of Ni a /Nu 
obtained in our study is larger than that concluded to be ~ 10% by Tsujimoto et al. (1995), as can be seen in table 
3. We think that this is the present situation of the predictability of numerical calculations, and that this dispersion 
must always be taken into account. For example, our results will bring the dispersion of both the average lifetime 
of SNe la's progenitors and history of the star-formation rate in our galaxy. Moreover, the critical mass, which 
represents the upper-mass limit of SNe II progenitor, may also be changed. 

A realistic s-process calculation will increase the amount of some less- produced elements in Ha95, as is shown 
by Prantzos et al. (1990). It would hold true for WW95. This suggests that the neutron-capture process during 
stellar evolution is also important for some elements. Since WW95 used a larger network for the 'post processing', 
they obtained neutron-rich elements. On the other hand, Ha95 used a small network during stellar evolution. As 
a result, neutron- rich nucleosynthesis seems to be insufficient. We also note that there is a possibility that the 
chemical composition of SNe la could be changed if the mass-accretion rate is changed (Nomoto et al. 1997b). An 
asymmetric explosion in SNe II may also help to solve these problems, since it would produce more 44 Ca and Ti than 
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the spherical calculations (Nagataki et al. 1997a). The initial conditions of the shock wave should also be considered 

(Nagataki et al. 1997b) because they affect the production of the iron-peak elements. 

Finally, let us to say that there remains still some physical uncertainties to be studied in order to clarify the 

uncertainty of supernova nucleosynthesis, which could be greater than the dispersion obtained in the present study. 

Therefore, we must trace the physical uncertainty in order to find the reason for the over/under production problems. 

Such efforts will cause the dispersion of Ni a /Nu to be smaller, and we should then be able to discuss the chemical 

evolution in the universe more precisely. 

We are grateful to an anonymous referee for useful comments. We also thank for Ken'ichi Nomoto for useful 
discussions. This research has been supported in part by a Grant-in- Aid for the Ccntcr-of-Exccllcncc (COE) Research 
(07CE2002) and for the Scientific Research Fund (05243103, 07640386, 3730) of the Ministry of Education, Science, 
Sports and Culture in Japan and by Japan Society for the Promotion of Science Postdoctoral Fellowships for Research 
Abroad. 
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Table 1. Nuclei used in the analysis. 



Species 


Case 1 


Case 2 


Case 3 


Case 4 


Case 5 


LMC 


SMC 


W7 


WDD2 


Ha95 


WW95 


lt)Q 


o 


o 


o 


o 


o 


o 


o 


1.43E-01 


6.93E-02 


1.80E-00 


1.14E-00 


1»Q 


o 


o 




o 


o 


o 


o 


8.25E-10 


4.62E-07 


4.61E-03 


4.04E-03 


90 H T 

20 Ne 


o 


o 


o 


o 


o 


o 


o 


2.02E-03 


9.13E-04 


2.12E-01 


1.65E-01 


91 H T 

21 Ne 


o 


o 




o 


o 


o 


o 


8.46E-06 


1.47E-06 


1.08E-03 


6.07E-04 


99 IV T 

22 Ne 


o 


o 




o 


o 


o 


o 


2.49E-03 


1.96E-06 


1.83E-02 


1.88E-02 


23 Na 


o 


o 




o 


o 




o 


6.32E-05 


1.30E-05 


6.51E-03 


5.13E-03 


Mg 


o 


o 


o 


o 


o 


o 


o 


8.50E-03 


4.76E-03 


8.83E-02 


3.83E-02 


9 "IV X 


o 


o 




o 


o 


o 


o 


4.05E-05 


2.39E-05 


1.44E-02 


9.15E-03 


9Rn .•■ 

26 Mg 


o 


o 




o 


o 


o 


o 


3.18E-05 


3.57E-05 


2.01E-02 


1.18E-02 


97 . i 

AI 


o 


o 




o 


o 






9.86E-04 


2.74E-04 


1.48E-02 


7.96E-03 


9S . 

28 Si 


o 


o 


o 


o 


o 




o 


1.50E-01 


2.71E-01 


1.05E-01 


1.14E-01 


9Q r-H . 

29 Si 


o 


o 




o 


o 




o 


8.61E-04 


3.87E-04 


8.99E-03 


3.61E-03 


30 Si 


o 


o 




o 


o 




o 


1.74E-03 


6.35E-04 


8.05E-03 


4.02E-03 


_ 

oi p 


o 


o 




o 


o 






4.18E-04 


1.80E-04 


1.21E-03 


1.14E-03 


^9 i-i 

32 S 


o 


o 


o 


o 


o 


o 


o 


8.41E-02 


1.65E-01 


3.84E-02 


5.34E-02 


33 g 


o 


o 




o 


o 


o 


o 


4.50E-04 


2.49E-04 


1.78E-04 


3.18E-04 


34 g 


o 


o 




o 


o 


o 


o 


1.90E-03 


2.50E-03 


2.62E-03 


4.33E-03 


35 CI 


o 






o 








1.34E-04 


9.83E-05 


1.01E-04 


3.80E-04 




o 






o 








3.98E-05 


3.36E-05 


1.88E-05 


7.59E-05 


3b Ar 


o 


o 




o 


o 


o 


o 


1.49E-02 


3.35E-02 


6.62E-03 


9.05E-03 


38 Ar 


o 


o 




o 


o 


o 


o 


1.06E-03 


1.45E-03 


1.37E-03 


2.22E-03 


39 K 


o 






o 








8.52E-05 


9.00E-05 


6.23E-05 


2.38E-04 


K 


o 






o 








7.44E-06 


7.12E-06 


5.07E-06 


5.31E-06 


40 Ca 


o 


o 




o 


o 


o 


o 


1.23E-02 


3.45E-02 


5.77E-03 


6.53E-03 


44 Ca 


o 






o 


o 


o 


o 


8.86E-06 


2.07E-05 


5.53E-05 


2.49E-05 


46 T j 


o 






o 








1.71E-05 


1.76E-05 


7.48E-06 


2.60E-05 


4 'Ti 


o 






o 








6.04E-07 


1.24E-06 


2.11E-06 


4.94E-06 


48 Ti 


o 






o 








2.03E-04 


8.53E-04 


1.16E-04 


3.17E-05 


4Q m- 

49 Ti 


o 






o 








1.69E-05 


6.71E-05 


5.98E-06 


4.10E-06 


50 Ti 


o 






o 








1.26E-05 


4.60E-04 


3.81E-10 


6.39E-06 


^O y-"1 

50 Cr 


o 


o 




o 


o 


o 


o 


2.71E-04 


5.24E-04 


4.64E-05 


8.16E-05 


^9 ^ 

52 Cr 


o 


o 




o 


o 


o 


o 


5.15E-03 


2.01E-02 


1.15E-03 


2.32E-04 


Or 


o 


o 




o 


o 


o 


o 


7.85E-04 


2.26E-03 


1.19E-04 


2.52E-05 


54 Cr 


o 


o 




o 


o 


o 


o 


1.90E-04 


2.03E-03 


2.33E-08 


1.54E-05 


55 Mn 


o 


o 




o 


o 


o 


o 


8.23E-03 


1.88E-02 


3.86E-04 


1.94E-04 


54 Fe 


o 


o 




o 


o 


o 


o 


1.04E-01 


7.08E-02 


3.62E-03 


6.52E-03 


56 Fe 


o 


o 


o 


o 


o 


o 


o 


6.13E-01 


6.15E-01 


8.44E-02 


1.66E-02 


57 Fe 


o 


o 




o 


o 


o 


o 


2.55E-02 


1.39E-02 


2.72E-03 


5.90E-04 


59 Co 


o 






o 








1.02E-03 


8.60E-04 


7.27E-05 


1.43E-04 


58 Ni 


o 


o 




o 


o 


o 


o 


1.28E-01 


3.34E-02 


3.63E-03 


7.93E-03 


60 Ni 


o 


o 




o 


o 


o 


o 


1.05E-02 


4.15E-03 


1.75E-03 


9.18E-04 


62 Ni 


o 


o 




o 


o 


o 


o 


2.66E-03 


1.36E-03 


5.09E-04 


2.94E-04 


63 Cu 


o 






o 








1.79E-06 


4.25E-05 


8.37E-07 


2.61E-05 


65 Cu 


o 






o 








6.83E-07 


1.59E-05 


4.07E-07 


7.47E-05 


64 Zn 


o 






o 








1.22E-05 


6.97E-07 


1.03E-05 


3.93E-05 


66 Zn 


o 






o 








2.12E-05 


3.18E-05 


8.61E-06 


1.23E-04 



O Nuclei used in the analysis of g(r). 
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Table 2. M n /M Ia for each model. 



Case W7+Ha95 W7+WW95 WDD2+Ha95 WDD2+WW95 



46 nuclei 1.85 1.23 1.79 1.19 

31 nuclei 1.86 1.24 1.80 1.20 

6 nuclei 2.33 1.53 2.07 1.36 

20 elements 1.85 1.23 1.79 1.19 

14 elements 1.86 1.24 1.80 1.20 



Table 3. The most proper value of r p . 



Case W7+Ha95 W7+WW95 WDD2+Ha95 WDD2+WW95 



46 nuclei 
31 nuclei 
6 nuclei 
20 elements 
14 elements 
N la /N u 



(2.8E-01) 

9.2E-02 

8.5E-02 
(2.9E-01) 

8.1E-02 
1.3E-01 (6.1E-01) 



1.1E-01 
9.2E-02 
1.0E-01 
1.1E-01 
9.6E-02 
1.0E-01-1.4E-01 



(1.4E-01) 

5.0E-02 

7.6E-02 
(1.8E-01) 

6.4E-02 
7.7E-02-(3.3E-01) 



6.1E-02 
5.4E-02 
9.2E-02 
7.7E-02 
6.3E-02 
5.6E-02-1.1E-01 



Parenthesis represents that g(r) is insensitive to r. In that case, r p is hard to be determined. 
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Tabic 4. The most proper value of r p /g(r p ) for SMC and LMC. 



Object 


W7+Ha95 


W7+WW95 


WDD2+Ha95 


WDD2+WW95 


LMC 


1.6E-01/5.2E-02 


1.6E-01/6.5E-02 


1.1E-01/3.6E-02 


9.6E-02/4.7E-02 


SMC 


1.9E-01/4.7E-02 


1.9E-01/4.1E-02 


1.5E-01/3.7E-02 


1.2E-01/2.9E-02 



10- 2 r ° = WDD2/W7 n 

10- 3 r o • = WW95/Ha95 1 

10-" E— , — 1 — , — , — , — 1 — , — , — , — 1 — ^_ 

30 40 60 

Mass Number 

Fig. 1.. 
Figure Captions 

Fig. 1. Comparison of chemical compositions. The open circles represent WDD2/W7 and the dots denote 
WW95/Ha95 

Fig. 2. g(r) for each model. The thick-solid, long-dashed, short-long-dashed, short- dashed, and dot lines correspond 
to Case 1, 2, 3, 4, and 5, respectively. 

Fig. 3. Comparison with solar values. The open circles represent Ha95+SNe la and the dots denote WW95+SNe 
la. Upper: W7 is used for SNe la. Lower: WDD2 is used. 
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Fig. 2.. 




Fig. 3. 



